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Uncertainties in Theoretical He I Emissivities: H II Regions, 
Primordial Abundance, and Cosmological Recombination 
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ABSTRACT 

A number of recent works in astronomy and cosmology have relied upon theoretical He I 
emissivities, but we know of no effort to quantify the uncertainties in the atomic data. We an- 
alyze and assign uncertainties to all relevant atomic data, perform Monte Carlo analyses, and 
report standard deviations in the line emissivities. We consider two sets of errors, which we 
call "optimistic" and "pessimistic." We also consider three different conditions, corresponding 
to prototypical Galactic and extragalactic H II regions and the epoch of cosmological recom- 
bination. In the extragalactic H II case, the errors we obtain are comparable to or larger than 
the errors in some recent Yp calculations, including those derived from CMB observations. 
We demonstrate a systematic effect on primordial abundance calculations; this effect cannot 
be reduced by observing a large number of objects. In the cosmological recombination case, 
the errors are comparable to many of the effects considered in recent calculations. 
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1 INTRODUCTION 

The need for accurate theoretical helium emission predictions is 
present in many aspects of astronomy and cosmology. These in- 
clude estimates of the primordial helium abundance, calculations 
of the cosmological recombination history, and, most recently, esti- 
mates of the time variation of the Higgs vacuum expectation value 
(GaBner, Lesch, & Arenhovel 2008). However, these works may 
involve underestimates of the uncertainties in He I emissivities. 

Standard Big Bang Nucleosynthesis (SBBN) yields are gen- 
erally a function of a single parameter, the primordial baryon-to- 
photon ratio, rj. An important test of SBBN is that, for a single 
value of this parameter, the predicted light nuclei abundances agree 
with the abundances deduced from observations. This is usually 
expressed via a concordance diagram (see Kirkman et al. 2003). 
CMB observations have tightly constrained r) (Spergel et al. 2007). 
However, the primordial He abundance, Yp, deduced from observa- 
tions sometimes disagrees with the SBBN yield. This is troubling 
because He is the most abundant of these nuclei. This is an area 
of active research (Fukugita & Kawasaki 2006; Peimbert, Luridi- 
ana, & Peimbert 2007; Izotov, Thuan, & Stasinska 2007, hereafter 
FK06, PLP07, and ITS07, respectively). The usual method for de- 
ducing Yp is to construct a relation between metallicity (typically 
taken as the oxygen and/or nitrogen abundances) and the helium 
abundance in selected extragalactic objects. The dY /dZ relation 



E-mail: rporter@pa.uky.edu 



is then extrapolated to zero metallicity (at which point all nuclei 
abundances should be equal to their primordial values). For a recent 
review of the subject, see Steigman (2007). Many authors have dis- 
cussed the errors involved in the dY /dZ relation analyses (Skillman 
et al. 1998; Peimbert et al. 2003; Porter et al. 2007). 

The study of the cosmological recombination spectrum is also 
currently an area of active research (Kholupenko, Ivanchik, & Var- 
shalovich 2007; Chluba & Sunyaev 2008; Wong, Moss, & Scott 
2008). The majority of effects considered in recent years modify 
the electron fraction by < 0.1% (see Switzer & Hirata 2008). 

In this Letter we seek to address the errors in theoretical He I 
emissivities. We have previously calculated He I emissivities with 
two different codes (Bauman et al. 2005; Porter et al. 2005, 2007). 
In the low-density limit Bauman et al. (2005) compares the results 
of these two calculations. We note that the differences listed there 
are not meant to be taken as errors. Here we analyze available in- 
formation and attempt to quantify the accuracy of all atomic data 
involved in our model atom (Porter et al. 2005). We then quan- 
tify the uncertainties in emissivities via two methods: 1) a rigorous 
propagation of error in the low-density limit; and 2) random pertur- 
bations in a Monte-Carlo analysis. We consider three sets of con- 
ditions typical of Galactic H II regions, extragalactic H II regions, 
and the epoch of cosmological recombination. Finally, we demon- 
strate that, while our analysis involves only statistical errors, the er- 
rors will yield systematic effects in primordial helium calculations. 
The most recent calculations of Yp (PLP07, ITS07) claim uncer- 
tainties < 1%, including estimates of the atomic data uncertainties. 
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Table 1. Assumed Uncertainties in Helium Atomic Data. 



Conditions 


Optimistic 


Pessimistic 


Rad. Recomb. Coefficients (Direct) 


n >=5 andL> 3 


0% 


0.1% 


n >= 5 andL <= 3 


0.01-0.7% 




n<5 


0.01-0.7% 


^4% 


El Transition Probabilities 


n„,ni < 10 and L < 7 


0.01% 


0.2% 


n„,ni < 10 and L >= 1 


0% 


0.01% 


n„ > 10, ni < 5 and L/ <= 2 


0.02% 


0.2% 


n„ > 10, L„ ^ 2 and L/ 5s 2 


0.6% 


4% 


n„ > 10, ni < 10 other 


1% 


7% 


«„,«; > 10 


10% 


10% 


Other Transition ProbabiUties 


2p3pi -Is'S 


1% 


5% 


2p3p2-l4'S 


1% 


1% 


2i^S - li'S (2v) 


10% 


30% 


2s^S-\s^S(Ml) 


1% 


20% 


2s^S- li'S 


1% 


5% 


all others 


1% 


1% 


Collisional Deexcitation 


rii, <= 5 and n/ <= 2 


10% 


30% 


A« = 


20% 


30% 


otherwise 


20% 


30% 



We demonstrate that atomic data uncertainties alone introduce a 
systematic error of similar size. 



2 ATOMIC DATA 

The physical processes discussed here are largely described in Os- 
terbrock & Ferland (2006). We group the relevant atomic data by 
physical process and present them in order of generally increas- 
ing uncertainty as follows: bound-bound radiative transitions, free- 
bound radiative transitions, and collisions. Additional sources of 
uncertainty are also considered in Section [2!4l below. Level ener- 
gies are by far the most accurate data involved, and we neglect these 
uncertainties entirely. 

Because some atomic data uncertainties will have a nearly lin- 
ear effect on corresponding emissivity uncertainties, it is critically 
important what uncertainties we assign to these data. In an attempt 
to mitigate this, we consider both pessimistic and optimistic cases, 
which are meant to respectively represent the minimum and maxi- 
mum uncertainty in each datum. In all cases we allow only one sig- 
nificant figure in our uncertainties, with optimistic cases rounded 
down and pessimistic cases rounded up. 



2.1 Transition probabilities 

Electric-dipole (El) transition probabilities between low-lying lev- 
els are extremely accurate. Drake (1996) calculated values for 
« ^ 10 and L ^ 7. We will refer to these levels as the "Drake set." He 
claims "essentially exact results for the entire singly-excited spec- 
trum of helium." For a small subset of these transitions, Argenti 
& Moccia (2008) compare their independent calculations with the 



corresponding Drake oscillator strengths; these agree to 5 or 6 sig- 
nificant figures. Argenti & Moccia also discuss discrepancies be- 
tween their results using velocity, length, and acceleration gauges 
and find that the acceleration gauge occasionally produces results 
that differ by one part in 10,000. We therefore conservatively adopt 
the still negligible uncertainty of 0.01% for the Drake transition 
probabilities. For the pessimistic case, we take 0.2%, an estimate 
of the higher order relativistic corrections given by Drake & Mor- 
ton (2007). 

Next we consider El transitions with upper level beyond the 
Drake set and lower level within the Drake set. These transitions 
have n„ > 11 or Z^(Y ^ 8 and tii ^ 10 or ^7. We use several 
different algorithms for these transitions. For ni ^ 5 and Z,/ ^ 2, 
we extrapolate Drake results as in Burgess & Seaton (1960). For 
transitions with both upper and lower L ^ 2, we use hydrogenic 
rates. For all other transitions with lower level within the Drake set, 
we use Drake's (1996) semi-classical algorithm for calculating ra- 
dial integrals from quantum defects. These algorithms agree with 
the tabulated Drake values to better than 0.05%, 4%, and 7%, re- 
spectively. On average, these respective algorithms agree to 0.02%, 
0.6%, and 1.0%. We adopt these maximum and average differences 
as the uncertainties of transition probabilities calculated with these 
algorithms in the pessimistic and optimistic cases, respectively. Re- 
maining El transitions have both upper and lower levels beyond the 
Drake set. Some of these rates (for example, with large n and small 
An) are expected to be significantly less accurate than the low-lying 
transitions calculated by Drake. We find, however, that because of 
the small contributions these levels have to the effective recom- 
bination coefficients of low-lying levels, the uncertainties in these 
transition probabilities are unlikely to contribute more than 0.01% 
to the uncertainties of low-lying emissivities. We assume 10% un- 
certainties for all of these transitions. 

Forbidden transitions are generally less accurate than 
their dipole-allowed counterparts. For the intercombination line 
we take the 1% band of theoretical values for 
the optimistic case and the roughly 5% experimental uncertainty 
for the pessimistic case (Dall et al. 2008). For the much slower 
2p-'P2 — li'S transition we take the transition probability from 
Lach & Pachucki (2001) and assume an uncertainty of 1% (al- 
though this transition should be completely negligible in all condi- 
tions). For the two-photon transition 2.s'S - 1^ 'S, we take 1% and 
5% for the respective optimistic and pessimistic cases. The former 
is roughly the dispersion in theoretical values (Drake 1996; Jacobs 
1971; Derevianko & Johnson 1997), while the latter is taken from 
the uncertainty in the experimental lifetime of 2s ' S given in Table 3 
of Derevianko & Johnson (1997). The two-photon El transition 
2j-'S - Is^S is assigned uncertainties 10% and 30%, based upon the 
standard deviation and spread in Table 4 of Derevianko & Johnson 
(1997). The corresponding Ml transition is assigned 1% and 20% 
based upon a discussion by Lach & Pachucki (2001). While more 
exotic transitions may be important in the cosmological recombina- 
tion context, their uncertainties should contribute negligibly to the 
uncertainties in our solutions. We assign 1% errors to these transi- 
tions. 

2.2 Free-bound radiative rates 

Photoionization cross-sections (used to calculate recombination co- 
efficients via the usual Milne relation) are calculated in a variety of 
ways which we will not repeat here. Uncertainties in the optimistic 
case are taken as follows. For « ^ 4, the threshold cross-section un- 
certainties are taken from the difference between the ab initio and 
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"extrapolated" results presented in Table 1 of Hummer & Storey 
(1998). For levels with larger n and L ^ 2, we simply use the same 
uncertainty that we assumed for 4^'^+'L. For levels with L > 2, 
we take the (already quite small) uncertainty assumed at L = 2. 
For simplicity, we assume in the optimistic case that the energy- 
dependence of cross-sections above threshold is exactly known. 
In this approximation, uncertainties in threshold cross-sections are 
equivalent to uncertainties in recombination coefficients, which we 
report in Table [T] 

In the pessimistic case, however, we apply uncertainties to 
recombination coefficients directly (rather than indirectly via the 
threshold photoionization cross-sections). For levels with « ^ 10, 
we set the uncertainty equal to the difference between the two codes 
discussed in Bauman et al. (2005). These differences reach as much 
as 4% at Te = 10^ K. We apply the results at n = 10 to all higher 
levels. Dielectronic recombination onto He+ forming He^ becomes 
important only at > 50,000 K, and in the cases we consider here 
it can be neglected entirely. 



2.3 Collision rates 

Collision rates are generally the least accurate data in the models 
(although they do not necessarily contribute the most to the total 
uncertainty in emissivities). Collisions between levels with An ^ 1 
are most efficiently driven by electron impact. Collisional deex- 
citation (and excitation) from levels with ;i„ ^ 5 and «/ ^ 2 are 
calculated using the close-coupling R-matrix results of Bray et al. 
(2000), and from clues therein we take 10% and 30% to be the 
optimistic and pessimistic uncertainties. For other transitions with 
An ^ 1 we employ various other algorithms and assume optimistic 
and pessimistic uncertainties of 20% and 30%, respectively. 

Angular-momentum changing (or Stark) collisions are most 
efficiently driven by slow-moving particles. In practice, these col- 
lisions are induced by protons. For non-degenerate transitions, de- 
fined for this purpose as transitions with An = and L ^ 2 and 
Al = ±1, we use Seaton (1962). In the energy-degenerate case we 
use the theory of Vrinceanu & Flannery (2001), which naturally 
treats Al > I transitions. While these "1-mixing" collisions are im- 
portant for driving highly excited states into statistical equilibrium, 
they are considerably less important for small n. We assume 20% 
and 30% for the optimistic and pessimistic uncertainties. 

Recent results suggest large uncertainties in Rydberg level col- 
lisional ionization rates. Nagesha & MacAdam (2003) found ex- 
perimental results that differed with respect to theoretical values by 
more than an order of magnitude. Vrinceanu (2005) has suggested 
additional physical processes are at work. Deutsch et al. (2006) ex- 
tended a separate theory of these collisions and also found results 
much less than the experimental results. Here we investigate the ef- 
fects of changing collisional ionization rates by large factors. Dis- 
abling collisional ionization entirely at = 10^ cm^-' causes line 
emissivities to increase by ~ 1% or less. Multiplying collisional 
ionization rates by 10 causes emissivities to decrease by ~ 1% or 
less. The three cases we will consider in Section |4] are each sig- 
nificantly less dense than this test case, so the effects of collisional 
ionization are even less important. We also note that collisional ion- 
ization should affect the Rydberg levels of He" and H" in similar 
ways. This means that any uncertainty in He I emissivities due to 
collisional ionization will be counteracted, at least to some degree, 
by uncertainty in H I emissivities when He I line fluxes are mea- 
sured relative to H I lines. We conclude that collisional ionization 
uncertainties are completely negligible. 



2.4 Other uncertainties 

The effect of the mixing of the singlet and triplet levels was con- 
sidered by Bauman et al. (2005). The effect was found to be neg- 
ligible for multiplet emissivities in the low-density limit. This is 
expected due to the principle of spectroscopic stability (Condon 
& Shortley 1991). That result may not hold at finite densities. To 
investigate this question, we use a mixing algorithm based on the 
method outlined by Drake (1996) and find effects that are negligi- 
ble in comparison with the assumed recombination coefficient er- 
rors. We note that Rubiiio-Martin, Chluba, & Sunyaev (2008) find 
interesting fine-structure absorption features in their y-resolved cal- 
culation of the cosmological He I recombination spectrum. 

The problems involved in modelling an (in principle) infinite 
set of levels with a finite system has been discussed at length by 
Porter et al. (2005) and Bauman et al. (2005). Because we need to 
run many models for the Monte Carlo analysis, we use a smaller, 
less computationally intensive model atom than was used in our 
previous works. We resolve all nLS terms with n ^ 40. While 
this smaller model yields emissivities that differ by nearly a per- 
cent from our larger model for some low-lying transitions, we per- 
formed a small Monte Carlo calculation with the larger model and 
found the dispersion in results comparable to the dispersion with 
the smaller atom. 



3 LOW-DENSITY, CASE B LIMIT 

In the collisionless, Case-B scenario, emission coefficients can be 
calculated without the inversion of a rate matrix. Rather, level pop- 
ulations are calculated by considering decay and cascade probabil- 
ities Pi,i and C„/ (defined following Robbins 1968). We use the un- 
certainties assigned in Table[T]and perform a rigorous propagation 
of error, neglecting any possible covariances not detailed above. 
The uncertainty in the emission coefficient 47lJ/nenif^+ is given by 

oj = hvui ^Jia^effP,if + {ap^,af^)^ (I) 

where a^^^ is the ejfective recombination coefficient (including in- 
direct recombinations from higher bound levels). 

As mentioned above. El transition probabilities, and therefore 
decay probabilities, are well-known for transitions between levels 
with low n. Thus equation[T]reduces to a linear relationship between 
Oj and C^i^'ff, and the uncertainties in dipole emission coefficients 
(in this collisionless, Case-B scenario) are entirely due to uncer- 
tainties in the effective recombination of the initial, upper level. 
This low-density analysis serves mostly as a check on our Monte 
Carlo results. Our low-density results are similar to the results of 
the 'extragalactic' Monte Carlo models discussed below. 



4 MONTE CARLO CALCULATIONS 

We define our errors as one-standard-deviation uncertainties. We 
disallow errors greater than three standard deviations from the nom- 
inal value because large deviations will inevitably produce unphys- 
ical trends with respect to quantum number. For the random number 
generator (RNG) we use an implementation of the very high peri- 
odicity Mersenne Twister (Matsumoto & Nishimura 1998), and we 
uniquely seed each calculation. We perform 1000 independent cal- 
culations for each model discussed below and for each set of error 
conditions in Table [T] 

Our three cases are defined as follows: 1) 'Galactic' - a case 
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Table 2. Standard deviations (%) of emissivities in our Monte Carlo analyis. 
Some trailing zeros are present for ease of comparison. 



Galactic Extragalactic 



X 


111 "+'L 




Opt. 


Pess. 


Opt. 


Pess. 


3965 


4p'P 


2i'S 


0.3 


0.8 


0.07 


0.2 


4471 




2p3p 


0.4 


1.0 


0.07 


0.4 


5876 


3d^D 


2p3p 


1.0 


4.0 


0.20 


0.8 


6678 


3d'D 


2p'P 


0.5 


2.0 


0.06 


0.3 


7065 


3s^S 


2p3p 


6.0 


20.0 


0.80 


3.0 


10830 


2p3p 


2i3s 


8.0 


30.0 


2.00 


9.0 


11970 




3p3p 


0.2 


0.9 


0.05 


0.4 


12530 


4p3p 


3*^8 


0.8 


2.0 


0.09 


0.4 


19540 


4p3p 


3rf3D 


0.8 


2.0 


0.09 


0.4 


20580 


2p'P 


2.V ' S 


5.0 


10.0 


0.30 


1.0 



B (Baker & Menzel 1938) calculation with Te = 10,000 K and 
He = 10"* cm"^ (similar to Models I and II of Porter et al. 2007); 2) 
'Extragalactic' - same as the Galactic case but with = 15,000 K 
and rie = 10^ cm^-'; and 3) 'Cosmological' - the He I recombina- 
tion epoch. The actual recombination history is known to be both 
non-time-steady and non-LTE (Seager, Sasselov, & Scott 2000). 
Here we consider the much less computationally intensive time- 
steady case. This model represents snapshots of recombining He+ 
in an extremely intense radiation field with cosmological parame- 
ters from Switzer & Hirata (2008a) and Sobolev optical depths. 



5 RESULTS AND DISCUSSION 

In Table (2] we present the standard deviation of the Monte Carlo 
emissivities for a number of emission lines. Note that one pair of 
lines is from the same upper level, 4p^P. We obtain the necessary 
result that these lines have identical uncertainties. 

In Figure [T] we investigate the effects of the present statistical 
analyses on Yp estimates. Each line represents a single full set of 
random perturbations (corresponding to a particular seeding of the 
RNG). For each set we calculate 7(A5876) for the range of tem- 
peratures shown on the x-axis. To first order all curves are linear 
offsets of the unperturbed line. This implies that statistical errors 
in He I emissivities would affect an ensemble of He abundances 
all in the same sense and to roughly the same degree, yielding a 
systematic offset in any dY /dZ relation. 

Figure [2] plots the SBBN helium yields and errors (solid and 
dashed lines, respectively. Buries et al. 2001) as a function of t]io 
(=1010t]). The blue boxes are the one standard-deviation results of 
a few recent Yp estimates. The higher and lower ITS07 boxes are us- 
ing the Porter et al. (2005) and Benjamin, Skillman, & Smits (1999) 
emissivities, respectively. The vertical yellow bar depicts the rjn) 
(and associated error) derived from 3-Year WMAP observations 
(Spergel et al. 2007). The red and green bars have heights equal 
to 2a (±1(7), where CT are present results for optimistic and pes- 
simistic extragalactic A5876 emissivities, respectively. The green 
bar could be reduced with a full set of high quality recombina- 
tion coefficients. It could be quite difficult to convincingly narrow 
the optimistic results represented by the red bar, thus prohibiting a 
dY / dZ-bsLsed derivation of 7] more accurate than those obtained by 
WMAP. 

PLP07 include a systematic error in Yp due to uncertainties in 
He I emissivities. Their value of 0.4% and their Yp error are consis- 
tent with the present uncertainties in extragalactic j(A5876). The 
FK06 error in Yp is also consistent with our results. Both the PLP07 



3 6 3 8 4 4.2 4.4 



Log (Tg/K) 

Figure 1. The emissivity of A5876 (multiplied by for easier viewing) 
for several random sets (dashed curves) of perturbed data versus electron 
temperature and the unperturbed result (bold solid curve). 




3 4 5 6 7 8 9 10 11 12 
Tijii (baryon to photon ratio) 

Figure 2. Primordial He mass fraction as a function of the baryon-to-photon 
ratio. Several recent results are indicated. See text for details. 



and FK06 results are consistent with SBBNAVMAP The ITS07 re- 
sults are inconsistent or marginally consistent with SBBNAVMAP, 
and their Yp errors are marginally consistent with our analysis. 

Correlated errors in pairs of He I emissivitites could be im- 
portant. We checked ratios of 114471, 5876, and 6678, in the op- 
timistic extragalactic calculations and found little correlation; the 
line-ratio errors are roughly the quadrature sum of the individual 
line errors. 

Finally, we consider the cosmological case. In Figure |3] the 
standard deviation in the electron fraction Xe = rie/rip is plotted 
against redshift. The pessimistic errors indicated by the green line 
reach 0.07%. The optimistic errors illustrated by the red line are 
roughly a factor of 3 smaller. Both are comparable to many of the 
effects considered by Switzer & Hirata (2008). A fuller non-time- 
steady analysis will be considered elsewhere (Porter et al., in prepa- 
ration). 

We thank Harry Nussbaumer for his helpful comments. 
RLP and GJF thank the NSF (AST 0607028) and NASA 
(NNG05GD81G) for support. 
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Figure 3. Fractional uncertainty in the electron fraction in the time-steady 
cosmological recombination model as a function of redshift. Red and green 
curves are the optimistic and pessimistic cases, respectively. 
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